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We study theoretically interspecies Cooper pairing in a fermionic system with SU(2) xSU(6) sym- 
metry. We show that, with suitable unitary transformations, the order parameter for the ground 
state can be reduced to only two non- vanishing complex components. The ground state has a large 
degeneracy. We find that while some Goldstone modes have linear dispersions, others are quadratic 
at low frequencies. We compare our results with the case of SU(N). 
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I. INTRODUCTION 

Higher symmetry groups such as SU(3) play an es- 
sential role in our understanding of elementary particle 
physics 1 . However, in typical condensed matter systems, 
the internal symmetries are much simpler. Superfiuid 
4 He has no spin, whereas 3 He and electrons in solids have 
only spin 1/2. The situation changes with the advances 
in cold-atomic gases. We have already seen many exam- 
ples of bosonic systems with finite spins (> 1)2. There 
has also been much attention in fermionic systems with 
more complex internal structure 3 -, and more recently, in 
systems where the symmetry is higher than the usual 
spin rotational symmetries. Examples include the hid- 
den symmetry in spin 3/2 fermions 4 , approximate SU(3) 
symmetry for 6 Li near special external magnetic field 
values^,— - — , and the enlarged symmetry for atoms with 
finite nuclear but no net electronic spin 14-16 . 

With no net electronic spin, the hyperfine spin of an 
atom comes entirely from its nucleus. Since the interac- 
tion between atoms mainly arises from their electronic 
clouds, the interatomic interaction is then independent 
of the total spins of the atoms involved. (This is in 
contrast with the more general situation where the scat- 
tering length between two atoms in general depends on 
their total spin. See re&2). A very interesting system 
of this class has been studied recently experimentally, 
namely a mixture of 171 Yb and 173 Ybil. m Yb and 
173 Yb have nuclear spins 1/2 and 5/2 respectively. For 
the 171 Yb and 173 Yb mixture under discussion, the intra- 
species interaction is weak (the s-wave scattering length 
is — 0.15nm between two 171 Yb atoms and 10.55nm be- 
tween two 173 Yb atoms). However, there is a rather 
large inter-species attraction, with the scattering length 
a = ai7i_i73 — 30.6nm 17 . The Kyoto group has al- 
ready been able to cool this mixture much below the 
degeneracy temperature, thus raising the interesting pos- 
sibility of interspecies Cooper pairing in this system. 

Motivated by this, we study a two species fermionic 
system with interspecies attractive interaction. The weak 
intraspecies interaction is expected to only slightly mod- 
ify the quantitative details of the system and will be ig- 
nored. We shall take a SU(2) internal symmetry for the 
first species, and SU(6) for the second one, though our 



findings are immediately generalizable to SU(2/ + 1) with 
general half-integer /'s. The system is expected to un- 
dergo interspecies Cooper pairing. (We shall only con- 
sider weak attractive interactions, and therefore ignore 
possibilities of bound states involving three or more par- 
ticles, c./i^ii 3 -). The general order parameter is thus a 
2 x (2/ + 1) complex valued matrix, since pairing can 
occur between any internal state of the first species and 
that of the second species. We determine the structure of 
the order parameter for the ground state, employing the 
mean-field approximation. We show that, with a suit- 
able choice of basis, the order parameter for the ground 
state can be reduced to only two non-vanishing complex 
components. The ground state for this system has also 
investigated before 1 ^ in the case where the interaction 
depends on the total spin F of the two interacting atoms 
and therefore the Cooper pairs have a definitive spin F. 
We shall mention briefly the relation of our results to 
theirs. 

The ground state is found to possess a large degener- 
acy. We thus proceed to find the Goldstone modes of 
this system. We show that there are 2 x (2/ + 1) such 
modes, and we shall evaluate their dispersion relations. 
We find 4 linear modes, and 2 x (2/ — 1) modes which 
have quadratic dispersion at very low frequencies, but be- 
coming also linear as slightly higher frequencies. We also 
determine the physical variables coupling to each of these 
modes, thus indicating how these modes can be excited 
experimentally. 

Currently, the fermi temperature Tp in the experiment 
o& is around 200nK, thus the product |fcjra| « 0.36 (here 
kp is the fermi wavevector). While this allows a weak- 
coupling consideration as in here, the transition temper- 
ature is unfortunately low (T c /Tp ~ e~ \ 2k F«\ ~ 0.01). 
However, we hope that eventually the superfiuid state can 
be reached (perhaps using optical Feshbach resonances to 
enhance the interaction 17 ) and the physics discussed here 
be studied. Also, we shall see that many of our physical 
results are more general, and thus would be applicable in 
case other more favorable related systems can be found. 

This paper is organized as follows. In section |H] we 
consider the ground state properties, and in section IIIII 
we discuss the collective modes. Our results have many 
similarities but some differences with the more studied 
case of SU(N) fermi superfluids. We therefore compare 



2 



our results with this case in Section ITVl We conclude in 
section [V] 



II. GROUND STATE 

In this section we discuss the ground state properties, 
assuming mean-field theory. This requires that the scat- 
tering length between the 171 Yb and 173 Yb atoms be 
small compared with the interparticle distances, and also 
that the transition temperature is small compared with 
the fermi temperatures. Nevertheless, in below we shall 
argue that many of our results are qualitatively correct 
beyond mean-field approximations. 

Let us denote the annihilation operators for the 171 Yb 
and 173 Yb atoms by , and eg respectively, where k 

is the wavevector and X = ±i, v = — /, ...,/ denotes the 
internal states. The Hamiltonian H has two parts. The 
kinetic energy is given by 



Here A is a 2 x (2/ + 1) matrix with elements Aa „. a, 
Pi can easily be evaluated but for here it is sufficient to 
know that all /3i > 0, and a is positive above some tran- 
sition temperature T c and negative below it. We would 
like to find the form of A which minimizes below T c . 
For this, we notice that since all fli > 0, A must be 
such that these higher order terms are minimized for any 
given Tr [AAt]. Let us denote D\ = [£„ |A AiI/ | 2 ] 1/2 - 

M where M is a Her- 



We obtain AAt = J hl3^1Il 



mitian matrix, Tr[M] = 0, and when expanded as M = 
Mi<ti+ M2U2 + M3<T3 using Pauli matrices ci 2 3, we have 
M 3 = (Df /2 - D 2 _ 1/2 )/2, M 1 - iM 2 = E.Ai^A*^. 
We thus have Tr [AAt] = £)2 + Jfi_ We also find 
easily Tr [(AAt) 2 ] = (D 2 1/2 +D 2 _ 1/2 ) 2 /2+Tr[M 2 ]. Thus, 

for given Tr [AAt] ; this fourth order term would be min- 
imized if we choose 
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/i, m is the mass of the atoms (We ig- 



nore the small mass difference between 171 Yb and 173 Yb 
here) For simplicity we shall also confine ourselves to the 
case where the chemical potentials /1 of the two species to 
be equal (and there arc no chemical potential differences 
among the different hyperfine sublevels)^. The interac- 
tion Hi n t is given by 



Hint =9j2 4 + ,A C l*L ,, c -fe'_ ,u a P + ,x 



(2) 



where k± = k ± q/2. This is the most general inter- 
species s-wave interaction obeying SU(2)xSU(6) symme- 
try, g can be eliminated in favor of the scattering length 
a m-i73 but we shall not need this explicit relation here. 

Within mean-field theory we can replace Hint by an 
effective interaction 



and 



D 1/2 = D-1/2 = D 



(7) 



Hence, if we regard ^1/2, v an d A-i/2,u each as a (un- 
normalized) wavefunction of a spin / particle, then eq 
© requires that these two wavefunctions are orthogonal, 
whereas eq ([7]) shows that they are of equal magnitude. 
A possible choice of A satisfying eq ([6]) is one where 
all elements Aa,^ vanish except Ai/2,1/2 and A_ 1/2,- 1/2- 
Eq {7J then requires that their magnitude to be equal. 
We shall give an alternate explanation of eq (j6]) and © 
below »2fl Using similar reasoning as above, we can actu- 
ally see that in fact alH > 2 terms are minimized by the 
conditions eq ([SJ) and (|7J|. The free energy then becomes 
= 2(aD 2 + /3 2 D 4 /2 + ...), the same as the usual BCS 
theory for a two-component system except an overall ex- 
tra factor of 2. 

Now we return to the microscopic theory. Defin- 
*A,„ via A x ,u = Dx^x, v thus E„I*a^| 2 = 



ing 



1, the pairing term in H, 



c f 

yL r 

— K,1S 
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can be written as 
and thus can be in- 



k wjriting, we are sii 



terpreted as pairing between at and the state c - 

' k , A 



-fe,A 



A 
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where has to satisfy the self-consistent equation 

Aa,, = (-5) J2 < H,x c -k, v > ( 4 ) 



Before solving this mean-field Hamiltonian we apply 
the Ginzburg-Landau (GL) theory. The GL free energy 
has the form 

n = aTr [AAt] + ^ ^i Tr [(AAt)'] ( 5 ) 

l>2 



, . (Here, for convenience of 

wjriting, we are simply calling the particles by their corre- 
sponding operators). Eq ^ implies that the most favor- 
are orthog- 



able state is such that & - 



and c - 



-fc,l/2 -fc,-l/2 

onal to each other. This is physical reasonable, as then 
at , and at , does not have to compete with each 

k,l/2 fe,-l/2 

already 



other to pair with the c- atoms. With c - 



fc,±l/2 

defined above, we can introduce c r for v ^ ±i to 

— k,u * 

make a complete set, therefore a unitary transformation 



between c - and c 1 - operators: 

— k,i> —k.y 



^k 
— k,v 



— k.v' 



(8) 
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where U v y = ^\= v y for v = ±|, and U v y for 
v 7^ ±i can be arbitrary so long as the matrix U is 
unitary. The kinetic energy can be re-written as K — 

be written as £a [^H.A^.A + 5 -k,\ a k,\) ~ ^ 
Thus the pairing term is of the normal BCS form ex- 
cept the sum over A. Thus the standard BCS calcu- 
lations can be immediately applied. The quasiparticle 

1/2 

energies are thus Et x = [£| + D%\ with D\ play- 
ing the role of the energy gap for the A species. For 
the expectation values, we have (for zero temperature, 
to which we confine ourselves for the rest of the pa- 
per) < a %x c_t x >= 



< 4.a°S,a >=< 4,a c M > 



±|. We also 



= n°(E kjX ) = i (l - etc for A 

have < a Z\C--t v >— whenever v ^ A. Expecta- 
tion values involving the c operators can be obtained 
by the inverse transformation. We get < a^ x c_t > 



A,;/ 



2E, 



2E, 



The self-consistent equation eq (|4]) 



becomes A\ tV = (-g) J2k 



ft 2E k 



and thus, after multiply- 



ing by v and sum over j/, either D\ = 0, or 



-*)E 
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cal ground state energy. The system is thus characterized 
by a large number of Goldstone modes, with the mode 
frequency uj vanishing as the wavevector q approaches 
zero. These modes are associated with the fluctuations 
of the order parameter components 5A\ jU (q) away from 
their equilibrium values. We shall employ the kinetic 
equation approach^! to evaluate their dispersion. This 
method is equivalent to the random phase approximation 
in diagrammatic approaches. Though we would employ 
the weak-coupling approximation, we shall argue that 
many of our results remain qualitatively valid for general 
interaction strengths, provided that the broken symme- 
tries for the ground state remain the same as that found 
within the weak-coupling approximation. For simplicity 
we shall restrict ourselves to zero temperature. 

To simplify our notation we shall drop the tildes on 
the c operators, or equivalently take a reference state 
where U is the identity matrix. Physical variables can be 
obtained easily by applying the unitary transformation 
U. We list the different types of modes in turn: 

(case 1): Modes corresponding to 8A\ M with v =^ ±i: 
For definiteness we consider 8Ai/ 25 / 2 (q). Besides the 
mean- field pairing terms in eq ([3]), we include 



Thus D\ obeys the usual BCS gap equation. From the 
form of the Hamiltonian, it is obvious that the most fa- 
vorable state would have both D 1 / 2 and D_i/ 2 finite, 
and by eq ([9]), both attain the usual BCS value and 



thus equal, consistent with eq (0. The particles ct 



ft . V 



for v 7^ ±~ are not involved in pairing. They maintain 
the normal state energies and there are thus (2/ — 1) 
remaining fermi surfaces, and < ct c^ v >= /(£fc), the 
fermi function. 

The above can be readily generalized to higher symme- 
tries. For example, for SU(4) x SU(2/ + 1) with / > 3/2, 
then the ground state has an order parameter which, in a 
suitable basis, can be reduced to pairing only between a\ 
and c\ with the same A. There are (2/+1) — 4 = (2/ — 3) 
fermi surfaces remaining normal. 



5H = Y,5A 1/2m {q)a\ 



fc+,1/2 -fc_,3/2 



(10) 



where k± = k± |. The hermitian conjugate (h.c.) term, 
involving 8A^ 2 3 ^ 2 and other 8A\ >V turns out to be de- 
coupled from the equations below. 8Ai/ 23 / 2 (q) has to 
obey the self-consistent equation 



SA 1 



/2,3/2 



(q) = (-a)Yl < a k + ,i 



/2 C -k. 



,3/2 



> 



(1) 



(11) 



III. COLLECTIVE MODES 



A. Dispersion relations: weak pairing limit 

The ground state therefore has a very high degeneracy. 
Any choice of the unitary transformation U gives identi- 



where the superscript W denotes the first order fluctu- 
ation contribution. Its equation of motion can be easily 



obtained using the hamiltonian H = Hk 
We get 



H 



eff 



SH. 



J 



1 dt < a k+,l/2 C -k^,3/2 > 



(1) 



< at 



fc+,1/2 -fc_, 3/2 



> (1) +A 1/2 , 1/2 <c t _ s+i/2 c_ £ ^ 3/2 



> 



(i) 



4 



-<5A iA3 / 2 (<t) (< 



(0) -<ar 1/9 4 , 

fc+,1/2 fc+,1/2 



> 



(0) x 



(12) 



Here the superscript stands for equilibrium expectation values. Thus we need also the equation of motion for 



^ -fc+,1/2 -fc-,3/2 ^ ■ 



^ <C -fc + ,l/2 C -fe-,3/2> (1) = -(^ + -^-)< C U.l 



, l/2 C -fc-,3/2 >(1) + A l/2,l/2 < 



a fc + ,l/2 C -fc_,3/2 > 



(1) 



-8*i/2, 3 /2® < cl s+ A/2 ai +1/2 



(0) 



(13) 



obtaining thus a closed set of equations for < i/2 c -k 3/2 anc ^ ^ c £ i/2 C -fc 3/2 - > ^- Fourier transform 
and solving for < , n c t 



6 i/2°-fc 3/2 > ^ 1 '* an d inserting into eq ([H]), we get 



w - £fc_ + ^fc+ 



w — 6c_ — £'fc + 




(14) 



where we have also used eq ^ to eliminate the coupling constant g. The term in the square bracket is a pair 
susceptibility. We can also directly obtain it from evaluating the response function of < i/2 c -k 3/2 

>W to SH. 

The second term is for adding to the ground state a pair of c and a particles with energies and Ek + ■ This process 
occurs only when the final states are available, hence the factor 5 ^1 + JjJ~) (1 — /(£fc-)) m the numerator. Similar 
interpretation applies to the first term which stands for annihilation of the pair. Note that the fermi factor restricts 
< for this term, hence Ek + — £&_ = £* + + | > 0. We also note that eq (fT4|) converges in the ultraviolet^ 
One can easily check that, if q = 0, eq (fT4")) is satisfied with lu = 0, showing that we indeed has a Goldstone 
mode. We next search for a solution to eq (fl"4"|) for small q and lu. For this, we add to eq (fT4|) the vanishing quantity 

f ~ TIT 1 2lF~^ — ^ 2£T — )' anc ^ re pl ace the dummy variable fe by k + , so that fc_ — > k. Expanding the 

V k + k + k + J 

resulting equations in q and lu, taking angular average, we obtain 







= Aiui ■ 



Bi 



B 2 f- 

m 



(15) 



where A\ 



-k 2E k 



£-ik (2m J 



2E- 



i-/fe) 

(Efc-Kk) 
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Efc* 
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fc 2E k 



(Efe-«fc) 2 



Si 



( Ek—£k E k J 



andS^fEs 



Cfc 1 )b 2 f 2 fc 2 
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Formally the term A 2 lu 2 is small and thus can be dropped in the lu — > limit. However, we shall see that in the 
weak-coupling regime, the coefficient A\ is small, and hence this term need to be kept in general beyond some small 
frequency regime lo* which we shall define later. 

Before we discuss this weak-pairing limit in which we are principally interested, it is instructive to evaluate first 
the dispersion in the strong pairing limit, where n < and 3> A (c./i 2 -^). Here A stands for the value of 

|^i/2,i/2 1 = I A_!/ 2 ,— 1/2 1 in equilibrium. In this limit, we get B\ = — |Ai, while B 2 is negligible. A 2 = Ai/(8\fi\) and 

2 

hence its contribution is always negligible when u -C \/j,\. Hence we obtain ui = 4—, the energy of a free particle of 
mass 2m (see also discussions below). 

Now we return to the weak-pairing limit, where ^ > A. The expression A\ is explicitly particle-hole asymmetric, 
i.e., if we approximate the density of states near the Fermi surface by a constant N(Q) and replace ~^ N(0) J d£, 

1/2 1 /2 

A would vanish. Hence we must use the more accurate expression — ¥ m ( 2 ^) — J*^ ^ M _|_ JL^ _ Dividing the 
region of integration to |£| < fi and |£| > /Lt, one can show that the later is smaller than the former in the /i> A limit. 



The first contribution to A\ can be re-written as 



ro(2m^) 1/2 rp, 



2tt 2 



1/2 



1/2' 



1 

2A 3 



The quantity in the square bracket can be Taylor expanded, and since ( 1 



m(2m)i) 1/2 



iS (^^a) in the ^ » A limit. 



(£2 + A 2)l/2 



A^ for $ > A, we find 



2tt 2 2A 

is even under particle-hole symmetry. Hence, the dominant contribution in the |A| <C /x limit can be approxi- 

m(2m)i) 1/2 rp. rlc 1 1 m(2m)i) 1/2 1 

W 1 J-fj, fl S 2B (B+C) 2 2^^ 2A 2 "' 



mated as A2 
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h' 2 

For B\, we substitute ^ = £, + Mi which generates respectively particle- hole asymmetric and symmetric contribu- 
tions (note the extra factor of £ in the definition of Bi), and the latter one is associated with a large coefficient /i. 

This term gives the dominant contribution to B\. We obtain B\ = — | m ( 2 ^) — ^-K" where if is the dimensionless 

integral A 2 ^ f -g^ + -> | for ^ > A. Thus S a > A. We find that B 2 < S x in the /i > A limit. 

(When using a constant density of states, it turns out that the particle-hole symmetric part of B 2 yields exactly zero, 
and hence we also need to include the energy dependence of the density of states. B 2 is of order (A/ fj,) 2 smaller than 
Bi.) We obtain thus finally the dispersion 

2 2 

= A\u) + — — B— (16) 
A m 

where A\ = ^hi-^, B = This is a quadratic equation and can be solved explicitly, but the main results can also 
be obtained by simply comparing terms in eq (fl~5j) . When |^ <g; A 2 , that is, when vfQ *C Vpq* = A(~ln-^), where 
vp = -\/2/i/m is the fermi velocity, we have 



lit 
3 A 



2m 



(17) 



The frequency is thus quadratic in q, but with a coef- 
ficient much larger than 1/m. Alternatively, we can also 
write 



In 2 

3 

fi A 



(18) 



where £o = ^a" is a measure of the zero temperature co- 
herence length. We note that the factor within the square 
bracket [ ] is much larger than 1. If q ^> q* however, the 
dispersion becomes linear, with 



1 BA 



v F 



(19) 



The transition between these two regime occurs at q « 
<f , where u w uj* = A^ln^. 

The existence of this quadratic mode at small q and uj 
can be understood from gauge symmetry arguments. The 
equations of motion such as eq (|12|) and (IT3l) must remain 
valid under arbitrary gauge transformations (we leave out 
the subscripts k to simplify the writing here) a a — > a\e i6x 
and c y — > c„e^". Hence, the equation of motion for 
5A l / 2Z / 2 must decouple from those for other dAx^ and 
also all 5A* x ^s. For example SAy 23 ^ 2 transforms dif- 
ferently from SAi/ 23 / 2 under a gauge transformation of 
c 3 /2, and this difference cannot be compensated by any 
combinations of Ai/2,1/2 and A_!/ 2 .--i/2- It also cannot 
couple to <5A_ 1 / 2i 3/ 2 : for example A^ 1 y 2 1 y 2 <5A_ 1 / 2i3 /2 
and A^y 2 1 ^ 2 <5A 1 /2,3/2 transform in the same way under 
transformations of ai/ 2 and a_i/2, but they do not do so 
under transformations of c±i/ 2 . We thus obtain a scalar 
(not the determinant of a matrix, c.f. (case 2) and (case 
3) below) equation which relates uj and q (such as eq 
(1141) ). The small uj and q expansion of this equation 



must therefore be of the form = Aiuj + £?— , where 
Ai does not vanish unless there is another symmetry 
(here particle- hole), and the expansion in q begins with 
q 2 since our system is spatially inversion symmetric ) 23 ' 24 . 
Therefore, though our specific formulas assumed weak- 
coupling, the quadratic dispersion at small q and uj is 
more general, provided that the broken symmetries re- 
main the same as those found within our weak-coupling 
theory. The frequency scale uj* unfortunately is small in 
the weak-coupling regime. However, we also note that 
uj* and q* increase with A//i, thus the frequency and 
momentum range where would have this quadratic mode 
would therer increase with the strength of the attractive 
interaction. 

As seen from eq (fT2|) and (fT3| . this mode couples to 
"spin-flip" of the c-particles. [This coupling is allowed 
because, e.g., in eq (TT2l . < a l / 2 c i / 2 > and A 1 / 21 / 2 < 
c i/2 c 3/2 ■ > transform in the same way under gauge trans- 
formations] . The order parameter mode under discussion 
therefore can be excited by Raman pulses which inter- 
convert v — 1/2 and v = 3/2 hyperfine sublevels. Bragg 
scattering experiments have already been performed in 
fermi gases^. Though that experiment does not involve 
" spin- flips" , it seems that a generalization of the method 
there can also observe spin-waves and hence the order 
parameter collective modes here. 

Similar discussion applies to all v 7^ ±1/2. There are 
thus 2 x (2/ — 1) such modes. The mode labeled by 
<5 A a v couples to the observable < c 1 -, c z >, a 

— fc+,A — K— i u 

generalized "spin-density". 

One can also consider 5A* X v with v 7^ ±1/2. These 
are just the complex conjugates of the modes discussed 
above, (with frequencies opposite sign 25 ) and are not new 
physical modes. 

{case 2): 8A\ tV {q) with v = A: We can proceed as 
above. This variable couples to SA* Xu (— q). The alge- 
bra is identical with the ordinary BCS case for pairing 
between two species, and hence we would not display 
the details of our calculations here. We just remind the 
readers that, due to the coupling between 8A\^(q) and 
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SA* X u (q), the dispersion is obtained by setting the deter- 
minant of a matrix to be zero^ (c.f., (case 1) above). In 
the weak-pairing limit we obtained the mode frequency of 
the Anderson-Bogoliubov mode 2 ^ ui = ^q. This mode 

can also be interpreted in a similar manner as in the two 
component case. (In the strong-pairing limit we obtain 
again the Bogoliubov mode for bound boson pair, similar 
to2£, but we shall not go into that here). These modes 
couple to the densities fluctuations < c 1 - c r % >^ 

^ -k +: \ -k-,X 

and < a r iS - >^. There are two such modes, 

one for each choice for A. 

(case 3) 5A\ t —\(q): we find that it couples with 
SA*_ xx (-q) [A*1 K _ X 6A X ,_ X and A X ^5A*_ XX transform 
under the same way under gauge transformations.]. The 
equation of motion is analogous to (case 2). Thus again 
we have linear modes with ui = -7=9. These modes cou- 



ple to the (spin) densities < - x c -u -a > ^ an< ^ 
< a t .at , > (1) . There are two such modes, again 

K +' A fe— ,— A 

one for each choice for A. 

The two modes in (case 2) can also be viewed as 
one "density" mode and one "longitudinal spin" mode 
(though the "density" mode involves only the total den- 
sities of the A = v = ±1/2 particles and the "longitu- 
dinal" mode involves the differences between them, and 
the v ^ ±1/2 ones are not involved). The longitudinal 
spin mode and the two "transverse" spin modes in (case 
3) give altogether three linear spin modes with oj — y=9- 

We had not discussed the spin density fluctuations cor- 
responding to < c' -. c_t „, >^ for u, 1/ + ±1/2. 

They do not couple to the order parameter fluctuations. 
Within our approximation they are simply independent 
particle-hole pairs as in the normal state. 



B. bosonic limit 



For strong attractive interactions, our formalism above may not apply due to the appearance of multi-particle bound 
states. However, to gain better understanding of some of our results above, it is instructive to consider this limit 
assuming we only have tightly bound pairs between a and c particles. We would like to illuminate on the counting of 
the collective modes and the existence of quadratic versus linear modes. Readers who find these points already clear 
are invited to skip this subsection. 

In this limit the system can be described by bosonic fields V'a.i/, corresponding to the bound state between a\ and 
c„. It is simple to construct a theory for a bosonic condensate of this system. The Hamiltonian H = Hk + Hi nt can 
be written as the sum of the kinetic part 



H K = 



E 

r \,v 



(20) 



where m& denotes the mass of the atom-pair and fib denotes the chemical potential, and the interaction, the most 
general form of which obeying SU(2)xSU(6) symmetry reads 



,A 3 ^A 2 ,A 4 + ^A 1 ,A4^A 2 ,A 3 ][<5 l / 1 ^3^2,f4 + $ui,VA&u 2 ,u^\ ltU ^\^ il/2 i>\ 3 ,v 3 ^A 4 ,i>4 ■ (21) 

For the ground state, we replace the operators i^x.v by c-numbers ^ \^ v . The resulting energy reads 

= -At 6 Tr[** t ] +||(Tr[** t ]) 2 +Tr[* 1 4' t ** t ]| (22) 



where SI/ is a 2 x (2/ + 1) matrix with entries ^\ >v . This 
energy can be minimized with exactly the same proce- 
dure as in sec HU We obtain J2 V ^i/2,u^*_ 1 / 2 v = 
and £„ |*i/ 2 ,„| 2 = £„ |*-i/ 2 ,,| 2 (cf. eq © and Q). 
One possible possible solution to these two equation is 
^1/2,1/2 = ^ -1/2,-1/2 = ^ but with all other compo- 
nents zero. Minimizing E, we get ^ = ^-\^\ 2 - Again we 



have many degenerate ground states. 

The collective modes are also in one-to-one correspon- 
dence with the ones we found above in the weak-pairing 
regime. They can be found by standard Bogoliubov 
transformation. We just state the results. The fluctu- 

2 

ation 5^\^ for v ^ ±1/2 has dispersion oj = cor- 
responding to free particles. There are 2 x (2/ — 1) of 
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2 

these modes. ( S^* X l , has u — — ^7 but simply corre- 
spond to removal of a boson and are not new modes). 
These are analogous to (case 1) discussed above. The 
four variables S^fx.x and 5^* x A with A = ±1/2 are cou- 
pled (the equation determining the relation between to 
and q is then again obtained by setting the determinant 
of a matrix to zero.). These yield two Goldstone modes, 
one with velocity lu = cq (tu = c s q) corresponding to in 
(out-of)-phase oscillations between the A = 1/2 and the 
A = -1/2 components. (5*i/ 2 ,i/ 2 = +(-)5*_i/2 -1/2 
and 5^X/ 2 i/ 2 = +( — )<^-i/9 -1/2)- We obtain c s = -^c 
and c 2 = /ib/irib- These are the analogous to the modes 
in (case 2). Finally, the four variables 6^x,-X an d S^* x _ A 
with A = ±1/2 yield two other Goldstone modes again 
with uj — c s q. These are the analog of (case 3). Wc 
have thus in total one "density" mode and three "spin" 
modes as discussed above. In the present case however 
the "density" modes and "spin" modes have different ve- 
locities due to interactions among bosons. The finding 
that the "density" and "spin" modes have the same ve- 
locities (as well as that they are given by Vp/y3) in the 
weak-coupling limit in the last subsection is the result of 
the approximation employed there, and is not expected 
to hold in general. 

IV. COMPARISON WITH SU(N) 

We now compare our results with BCS pairing in 
SU(N) models , where there are N species of fermions 
&i,...,N with interspecies interaction which is invariant 
under SU(N) transformations. We have seen that, for 
our SU(2) x SU(6) system with interspecies pairing, for 
the ground state, the pairing order parameter reduces 
to, with a suitable choice of basis, one where each o- 
species pairs only with one unique c-species. We have 
seen that this is a consequence of eq $3§. The results for 
SU(N) are similar— £ but there is one important differ- 
ence. For SU(N), the order parameter A is an antisym- 
metric N x N matrix (there is no such restriction for our 
SU(2)x SU(6) case). Under a unitary transformation U, 
it transforms as A — » UAU* where the superscript t 
denotes the transpose. Hence, with a suitable U, A can 
always be transformed^ to one where all entries A Xl \ 2 
vanish except 

A12 = -A21 A 34 = -A43 ... , (23) 

that is, 

A13 = A 14 = ... = A 23 = ... = (24) 

so that 1 only pairs with 2, 3 only pairs with 4 eto^i. 
However, we emphasize that the origin of eqs (|23[) and 
(|24|) is very different from the SU(2) x SU(6) case. The 
possibility of writing A in the form of Eqs (|23l) and 
(IM|) is purely a consequence of the fact that A is 
antisymmetric^, which in turn only requires fermionic 



anticommutation relations, and therefore holds for arbi- 
trary states (including excited states) for the system. For 
our SU(2) x 5(7(6) system, eq ^ needs not hold other 
than the ground state. 

Now we compare the collective modes. Consider first 
N = 3. In the gauge where the only finite component 
is Ai 2 , it can easily be shown that 5A\3 obeys the same 
equations of motion as eq (|T2"j) and (flU)) with ai/ 2 — > 
oij c-i/2 a 2, c 3 / 2 — > 03. Hence the dispersion found 
there also applies. That <5Ai3 decouples from other order 
parameter modes can also be seen by gauge invariance 
arguments as in sec IIII1 It can be easily seen that 8A 2 z 
yields yet one more mode with the same dispersion as 
SAi 3 . On the other hand, <5Ai 2 and 8A\ 2 are coupled 
together, and they generate one linear mode as in case 2 
in sec IIIIl Hence in total, we have one linear mode, and 
two modes with quadratic dispersions for small uj and q. 

Our results for the number of linear and quadratic 
modes agree with Ref^. (note that they have an alterna- 
tive but related argument for the existence of quadratic 
modes). However, the precise dispersion of the quadratic 
mode is different. They obtain a dispersion same as a 
free particle of mass m. Our eq (|14[) is the same as 
their (Al) at zero temperature (except k — > k — |, see 
however—), so we believe that they have made a subse- 
quent algebraic error (the contributions from our Bi and 
B 2 terms seem to be missing). That their result is un- 
reasonable was also pointed out in RefJ^ (see also 2 ^) . If 
particle-hole asymmetric terms are ignored, we saw that 
the quadratic modes become linear with velocity vf/V^, 
in agreement with Ref J^, resulting in therefore three lin- 
ear modes in total. RefJ^ however did not provide the 
full dispersion relation in the presence of particle-hole 
asymmetric contributions. Also, we do not understand 
their claim that, when particle-hole asymmetric terms 
are included as we have done here, "two of the massless 
modes split into a massless mode with quadratic disper- 
sion and a massive one". We have seen that SAis and 
<5A 2 3 yield two independent quadratic modes. These two 
modes would not couple to each other, nor to other or- 
der parameter modes by gauge symmetry, and thus each 
must yield one quadratic mode. Since this is due to gauge 
symmetry, this result should hold even beyond the weak- 
coupling approximation we have employed. Lastly, Rcf. 7 
has counted the number of modes differently from here 
(and therefore also Ref. -). They regard <5A^ 3 and 5A 23 
as giving two other modes in addition to c)Ai3 and <5A 2 3, 
and they therefore counted in total five Goldstone modes. 
As remarked before, i5A^ 3 and <5Ai3 just correspond to 
annihilation and creation of the same excitation, so it 
seems natural not to count them separately (see sub- 
section IIIIBj) . Also, Refji has performed a numerical 
calculation of the dispersion for the mode <5Ai3. Their 
numerical result seems only to show a linear dispersion, 
likely due to the small A//J, values chosen there. 

The above discussions can readily generalized to larger 
N. For example, for SU(5), there are four quadratic 
modes. With choice of order parameter as in eq (|2"3"|) and 
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(124|) . they are 6A A ,5 for A = 1,2,3,4. No such modes 
are expected for N being even. For SU(4), there are six 
linear but no quadratic modes. 

V. CONCLUSION 

In this paper, we have considered some superfiuid prop- 
erties of an SU(2)xSU(6) system with interspecies pair- 
ing, motivated by the system studied experimentally in 
re&£. We considered both the ground state and collective 
excitations. Some properties are dramatically different 
from two-component systems. There are in particular 



collective mode excitations with quadratic dispersions at 
low frequencies. Many of our results found are generally 
applicable to systems with interspecies pairing with high 
symmetries. 
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5A ± 
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where the 
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!-/(«*_) /(f*_) 



-k 4 
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